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ABSTRACT 

We investigate the possibility of using globular clusters as targets for microlensing 
searches. Such searches will be challenging and require more powerful telescopes than 
now employed, but are feasible in the future. Although expected event rates are 
low, we show that the wide variety of lines of sight to globular clusters greatly en- 
hances the ability to distinguish between halo models using microlensing observations 
as compared to LMC/SMC observations alone. 
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1 INTRODUCTION 



The over 15 microlensing events observed towards the LMC 
represent a triumph of experimenta l perseverance in the face 
of seemingly insurmountable odds ( Axelrod 1997). However, 
microlensing has perhaps opened up more questions than it 
has resolved. What are the lenses? Where are they? What 
fraction of the halo do they constitute? Are they really part 
of the extended halo? Even if we could unequivocally state 
that the answer to the last question is "yes" , there are still a 
number of unanswered questions. Present observations can 
tell us very little about one of the most important issue: the 
structure of the baryonic halo. What is the shape of the dark 
baryonic halo? What is the core radius? Is the halo really a 
1/r 2 halo? Indeed, rather than being elucidated, these issues 
of structure cloud the interpretation of the events observed 
(Gates et al. 199(:). Studies show that this is inevitable even 



with a very large number of events ( Markovic fc Sommer- 



Larsen 1997). The difficulty in obtaining information on the 
halo structure stems from two sources. First, the basic mi- 
crolensing degeneracies prohibit obtaining direct distance in- 
formation along the line of sight. Second, even if such infor- 
mation could be obtained, inverting one dimensional data to 
obtain a three dimensional distribution is risky and prone 
to error. 

Two strategies can be employed to address these prob- 
lems: usin g parallax m easurements to break the basic de- 
generacy (Gould 1992) and using multiple lines of sight to 
extract structure information directly. A parallax satellite 
(or better a pair) would allow us to determine the distances 
to individual lenses and hence have information about the 
distribution of mass along the line of sight to the LMC. 
Given the expense of satellites, and the difficulties in trans- 
lating this to three dimensional information, it is reasonable 
to more fully explore the second possibility. This paper will 



look at the benefits of multiple lines of sight. They allow 
a more direct mapping of the variation of density with po- 
sition and thus a determination of such quantities as halo 
flattening. Ideally the two methods would complement each 
other allowing a full mapping of the halo with greatly less- 
ened theoretical uncertainty. 

A variety of possibilities come to mind when alternate 
lines of sight are contemplated: 

i) Andromeda (M31) 

ii) dSph 

iii) Globular clusters 

iv) bulge 

v) spiral arms 

Of these, the first provides only a single line of sight through 
our halo and of course will be complicated (form our point of 
view) by self lensing. The second possibility is more attrac- 
tive. Dwarf spheroidals have a large number of stars in an en- 
vironment where crowding is not expected to be a problem. 
Unfortunately they are both distant and few. Further, some 
are not even as luminous as large globular clusters. Discard- 
ing the most problematic, only a handful remain as likely 
candidates, reducing their value as probes of varied lines 
of sight. Finally, dwarf spheroidals are not well-understood. 
Their observed velocity dispersions seem to either require 
a substantial dark matter halo for each dSph or else indi- 
cate that they are extended structures fo rmed from the tidal 
disruption of a larger system (see e.g. (Gallagher & Wyse 



1994)). In either case, self- lensing may be an important con- 



sideration. 

The last two directions are expected (and observed in 
the case of the bul ge) to be dominated by lens ing from non- 
halo populations (Alcock et al. 1997a 1997a). Information 



about halo properties may be difficult to disentangle from 
the tail of the disk/bulge lensing distribution. Globulars are 
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well stu died systems, do not have halos of their own (Moore 



199f), i ,nd are small enough so that self lensing is negligible. 
Further, they are abundant, relatively close (but not too 
close) and contain a respectable number of stars. 

The rest of this paper is organized in four sections. 
The first discusses the selection of appropriate clusters and 
touches briefly on the expected number of stars observable. 
The next two sections set up characteristic galaxy models 
we will explore, discuss how effective cluster microlensing is 
in distinguishing between them and compare this to obser- 
vations of the SMC. Finally, we conclude with a summary 
of the expected benefits of cluster microlensing and what it 
will take to make it happen. 



2 SELECTION OF GLOBULAR CLUSTERS 

Observing globular clusters for microlensing will be dim- 
cult. The number of stars is limited both by the low masses 
of clusters and the highly crowded conditions. We discuss 
these challenges in more detail later, but for now it should 
be noted that the ability to observe lensing events in a glob- 
ular cluster is only half the story. The other half is whether 
such observations are worth the effort: will the information 
from observations along multiple lines of sight yield enough 
additional information as to the structure of the halo? Nev- 
ertheless, it is clear that not all globular clusters can be used 
and that the best candidates must be selected. 

We examined the globular cluster catalog made avail- 
a ble by William E. Harris (McMaster Univ ersity) at 
" tittp:/ / www.physics.mcmaster.ca/Globular.htmj " . To be 
useful for a microlensing search, a globular cluster must sat- 
isfy a number of criteria. First, it must have a reasonable 
number of potentially observable stars. Second, the optical 
depth must be relatively large, and third the microlensing 
rate due to halo objects must not be swamped by events 
due to disk and bulge lenses. We examine each of these con- 
straints in turn. 

Ideally, we would have simulated a microlensing search 
to each of the w 150 known globular clusters, creating images 
based on realistic cluster models (core radius, tidal radius, 
luminosity function, etc.) and observational parameters (see- 
ing, aperture etc). These simulated images could then be 
analysed with a photometric pipeline to determine the ef- 
fective number of stars observable for each cluster. Such a 
procedure would be very complex and time consuming. Fur- 
thermore it would be strongly dependent on the telescopic 
parameters and the observing strategies used. We adopted a 
very different procedure. Clusters were cut on distance mod- 
ulus and total luminosity. Estimates based on standard glob- 
ular cluster luminosity functions indicate that a My = — 7.5 
cluster should have ~ 200,000 stars above My = 10. We as- 
sume such stars can be observed out to a distance modulus, 
m_D, of 17.5. We discuss the challenges of observing stars 200 
times fainter than present microlensing searches in the our 
final section. We thus consider only clusters with both dis- 
tance modulus, mu <17.5 and My < —7.5. Many stars will 
be unobservable due to crowding. On the other hand, the 
typical cluster we select will be both brighter and closer, in- 
creasing the number of stars above the limiting magnitude. 
As a reasonable compromise we adopt 100,000 stars for each 



[hb] 

Table 1. Globular clusters that passed distance, brightness and 
latitude cuts. Distances from the Sun are in kpc. Visual lumi- 
nosity (My) and distance modulus (m^) are in magnitudes. The 
optical depth is for a standard halo and is in unit of 10 — 8 . 
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My 
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15.2 
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-7.68 
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6.4 
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1851 


11.7 


244.5 


-35.00 
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3.6 
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8.9 


282.2 


-11.30 
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11.3 
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24.6 


342.2 
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10.0 
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-7.90 
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15.8 
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-8.35 
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21.0 
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6712 


25.4 


5.6 
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21.5 
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15.27 


4.5 
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11.1 


53.4 


-35.80 


-8.90 


15.37 


5.9 



cluster. We discuss how our results change as this number 
is increased or decreased. 

We want our candidate clusters to be worth observing: 
Clusters with a very low optical depth in all models will give 
us little leverage for distinguishing between models. Since we 
do not know a priori what models might be interesting we 
used a halo independent measure, the heliocentric distance. 
Clusters close to the observer will have narrow, short mi- 
crolensing tubes and hence low optical depths. Conversely, 
distant clusters are much more likely to have high optical 
depths. Further, local clusters are much less useful as a probe 
of the global structure of the halo even if they have a high 
optical depth. We thus adopted a cut on the distance to the 
globular clusters, R s > 8 kpc. Clusters closer than this had 
negligible lensing. 

Finally, we wished to concentrate on the halo. Thus 
clusters for which non-halo lensing would dominate had 
to be removed. One could calculate the expected ratio of 
disk/bulge lensing versus halo lensing, but this procedure 
is fraught with the uncertainties of disk and bulge models. 
We took a simpler approach, eliminating those clusters lying 
within 10° of the galactic plane. 

Twenty clusters survived these cuts with an average 
luminosity My « —8.5 and distance modulus mn ~ 16. 
Optical depths for a standard halo model normalized to 
tlmc = 2.5 x 10~ 7 ranged from 0.3 to 4.6 xl0~ 7 . The clus- 
ters are listed in Table 1. 



3 MODEL DISCRIMINATION 

We know little about the detailed structure of dark halos. 
Rotation curves probe the radial profiles of galaxies, but it is 
difficult to determine the relative importance of the dark and 
luminous material. Further, rotation curves can only provide 
information in one direction. Flaring of HI gas layers in disks 
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Table 2. Characteristic models used. The core radii are given in 
kpc. 



Model 


Core 


Index 


Axis Ratio 




A 


5.0 


2.0 


1.0 


Standard halo model 


B 


20.0 


2.0 


1.0 


large core model 


C 


5.0 


2.0 


0.(5 


flattened model 


D 


1.0 


3.0 


1.0 


non-halo model 


E 








debris model 



would seem to be a good probe of the shape of the po tential 



whe n us ed in conjunction with the rotation curves (Oiling 
seem to give anomalous results that are difficult 
to reconcile with stability concerns or other tracers of the 
potentials. Polar ring galaxies provide a unique opportunity 
to probe the shapes of dark matte r halos in that ther e are 
orthogonal probes of the potential (Sackett et al. 1994), but 
the details of polar rings are uncertain and these galaxies are 
almost certainly deeply disturbed systems. Weak lensing has 
not yet borne fruit while simulations are not clear and could 
be lacking resolution and/or physics. 

Expectations for a baryonic halo are even more uncer- 
tain. A very rough idea of the uncertainties involved can be 
gained by adopting the following parametrized model for the 
halo: 



P = Po 



afi+ RP + {z/qY ' 



(1) 



Such a model has 4 unknowns corresponding to the mass, 
flattening, core radius and radial fall-off of the halo. Clearly 
the parameter space represented by such models is very 
large. Since the amount of data from cluster lensing will 
be small measuring precise values for the parameters is not 
expected. Instead, we will concern ourselves with how well 
cluster lensing can discriminate between classes of models. 
We pick representative models illustrating various portions 
of the parameter space. 

Our models are summarized in Table 2. The first model 
is the standard halo model of Griest (1991), with an isotropic 
velocity dispersion a ~ 156km/s. The second is a similar 
model, but with a much larger core radius. Third is a flat- 
tened halo as suggested by polar ring and HI flaring studies, 
with the velocity dispersion reduced in one dimension ac- 
cording to the tensor virial theorem. Fourth is a more steeply 
declining halo based on the stellar halo, with a velocity dis- 
persion given by a ~ 127km/s. This model is supplemented 
by a heavy thick disk which does about 30% of the lensing. 
As the lenses are not in the dark halo proper, we call this the 
non-halo model. Finally, fifth is a model with no MACHOs 
in the halo at all. The LMC lensing is provided by a tidal 
tail or a disrupted dwarf galaxy in the line of sight as sug- 
gested by Zhao (1997). All models are constructed to have 
tlmc = 2.0 x 10 -7 and Y LM c = 4.0 x 10 _14 s _1 . Microlens- 
ing observations towards the LMC alone will be incapable of 
distinguishing between these models since the distribution 
of lens masses is unknown. 

To determine how well cluster lensing can distinguish 
between models, we proceed in the following manner. First, 
we pick a "true" model S from the set described above. After 
computing the microlensing rate towards each target under 
this model we generate a realization of S. That is to say we 
pick an observed number of events for each target based on 



the rate and ten years of observation, using Poisson statis- 
tics. For each realization we compute the likelihood for each 
of our models. The most likely model will not necessarily be 
the "true" model. After generating many (100,000) realiza- 
tions we can build up the probability P(T) that the model S 
will be identified by a microlensing experiment as the model 
T. Repeating but with a different "true" model we build 
the full matrix P(S,T) which gives the chance that the real 
galaxy S will appear to the observers as T. If P(S,T) is di- 
agonal then we know that the observations are doing a good 
jobs of distinguishing models. If on the other hand, it is far 
from diagonal, then observations are doing a poor job and 
misidentifications are frequent: the data cannot even distin- 
guish between, for example, a flattened halo and a spherical 
halo. 

We present results from a number of scenarios: 
LMC/SMC lensing alone, LMC and globular cluster lensing, 
and LMC/SMC/cluster/dSph lensing. For the dSph lensing, 
we looked at lensing information that could be gained from 
looking at Fornax, Sculptor, Sextans and Ursa Minor. The 
other dSph's were either too far, too faint, or some combina- 
tion of the two to be useful. We take the effective number of 
stars in the SMC to be 2.0 x 10 6 . We also show how these the 
effectiveness of globular cluster lensing varies with number 
of observable stars. 



4 RESULTS 

It has been suggested that SMC lensing can break the degen- 
eracy between flattened and spherical halos (eg. Sackett & 
Gould 1993). Unfortunately, for realistic event rates the time 
required is large. As can be seen in Table 3, for a ten year 
experiment LMC/SMC observations can correctly select the 
flattened model only ~60% of the time. Examination of the 
likelihoods directly shows that even when the correct model 
is selected, the margin of probability is small with spherical 
models and large core radii models only slightly less pre- 
ferred. When the "true" model has a large core SMC lensing 
performs abysmally. Indeed, most of the time another model 
is preferred. As is true of the LMC, the line of sight towards 
the SMC is most sensitive at large galactocentric distances, 
where the effect of a larger core radius is small. As expected, 
SMC lensing does better for non-halo models. Here the di- 
rection of the SMC line of sight closer to the Galactic center 
allows it to penetrate into the dense inner regions. However, 
there is still a reasonable chance (30%) that an unexpect- 
edly high amount of SMC lensing mimics one of the other 
models. Further, examination of the likelihoods again shows 
that the probabilities are not sharply peaked on the non- 
halo models. Thus our confidence in this determination will 
not be high. The debris models are always correctly iden- 
tified. LMC/SMC lensing in general does not fare well in 
distinguishing between models because it only probes two 
lines of sight. This is particularly worrisome if s elf lensing 



is important as has been suggested for the SMC ( Palanque- 



Delabrouille 1997). Even a small amount of unknown self 
lensing is capable of throwing the model determinations into 
grave doubt. Furthermore, if the detailed structure of the 
baryonic halo is lumpy, as suggested in some models, the 
determinations become even more uncertain. Microlensing 
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Table 3. Discrimination matrix for five scenarios. See text and 
Tabic 2 for explanation of the models used. 

SMC only (2.0c6 stars) 
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E 




A 0.29 
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0.69 
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A 0.63 


0.04 


0.33 


0.00 


0.00 




B 0.03 


0.97 


0.00 


0.00 


0.00 




C 0.35 


0.01 


0.64 


0.00 


0.00 




D 0.00 


0.00 


0.00 


1.00 


0.00 




E 0.00 


0.00 


0.00 


0.00 
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SMC, G. Clusters and dSph 





A 


B 


c 


D 


E 


A 
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towards the SMC appears to be a poor tool for elucidating 
halo structure. 

Globular cluster microlensing can avoid some of these 
problems. First of all, quick calculations show that r se i f -C 
1CP 8 for all globular clusters. Thus, self lensing will be neg- 
ligible. Second, since many more lines of sight are involved, 
lumpiness in the detailed distribution of the mass will be 
largely smoothed over. The figures for LMC/globular clus- 
ter lensing shown in Table 3 should therefore be more robust 
than the corresponding figures for LMC/SMC lensing. 

A quick examination shows that globular cluster lensing 
is much better at distinguishing between models with differ- 
ent core radii. This is quite understandable: globular clusters 
probe precisely the region where differences in the core are 
most apparent. Similarly, the selection of the the non-halo 
models certain. Models with a large power-law index predict 



Figure 1. Distribution of best-fit power law indices for a true 
model of n=2.5 and LMC/SMC/cluster/dSph lensing. Clusters 
are assumed to have 100,000 stars monitorable. 



much higher lensing towards the globular clusters making 
the distinction very easy and certain. LMC/globular clus- 
ter observations do not fare better than LMC/SMC obser- 
vations at distinguishing flattened models. This is because 
only a few clusters are at high galactic latitude. As before, 
the debris models are firmly ruled out. We also examined 
a scenario where the SMC, globular clusters and four dSph 
galaxies were all observed. For these purposes we assumed 
that it was possible to monitor 500,000 stars in each dSph. 
Results are shown in the third part of Table 3. The improve- 
ment in the determination of the flattening is considerable 
due to the relatively high latitudes of the dSph's. 

Since the core radius and the power-law index are so 
easily distinguished by the full set of targets, we investigated 
the question of how well quantitative information could be 
obtained. We picked a standard model with n=2.5. Real- 
izations of this model were then compared to models where 
the optical depth towards the LMC was held fixed but the 
power law index was allowed to float. Figure 1 shows the dis- 
tribution of power law indices inferred. The peak is clearly 
around the "true" value with a small (~ 0.2) variance. The 
concentration of the globular clusters towards the Galactic 
Center allows very accurate determination of this parame- 
ter. A similar test was performed for the core radius, this 
time with a core of 12 kpc. Figure 2 shows the distribution 
of inferred values. Note that the distribution is not quite as 
tight (er ?s 5 kpc), but still allows fairly accurate determina- 
tion of the core radius. 

Finally, we examined how much our results depended 
on the number of stars observable per cluster. Table 4 shows 
the results of LMC/globular cluster lensing for 50,000 and 
200,000 stars per cluster. We see that the results are not 
strongly sensitive to increases in the number of stars ob- 
served, although marginal improvement is obtained. When 
the number of stars is decreased, however, the returns on 
cluster lensing are lower. Nevertheless, clusters can still dis- 
criminate between models better than LMC/SMC lensing 
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Core Radius (kpc) 

Figure 2. Distribution of best-fit core radii for a true model 
of a=12 kpc and LMC/SMC/cluster/dSph lensing. Clusters are 
assumed to have 100,000 stars monitorable. 



especially when the possible problems with LMC/SMC lens- 
ing are considered. Globular cluster microlensing is a worth- 
while investment, even at lower numbers of effective stars. 



5 CC NCLUSIONS 



Globular cluster microlensing searches are much better for 
determining important halo parameters than LMC/SMC 
measurements alone. Although few events are involved, the 
diverse lines of sight allow far greater leverage on halo struc- 
ture parameters. In particular, the power-law index and the 
core radius are readily accessible to cluster data. 

The importance of diverse lines of sight has already been 
demonstrated: reports of a microlensing event seen towards 
t he SMC by both the MACHO and EROS collabora tions 
( |Palanque-Delabrouille 1997| ; [Alcock et al. 1997c 1997c| ) sug- 
gest that the debris models for LMC lensing are not viable. 
Caution must be taken however, as the self lensing due to the 
SMC could cloud this interpretation. The small size and low 
self-lensing of the globular clusters is a definite advantage in 
this respect. 

The faintness of the typical globular cluster stars 
(Mv ~ 26), and the highly crowded conditions (>100 
stars/arcsec 2 at the core) both make microlensing obser- 
vations difficult. This suggests that the best bet for clus- 
ter microlensing studies would be a large (« 8m) telescope 
with excellent seeing (< 0.5"), possibly achieved with adap- 
tive optics. Telescopes of this type, such as the SUBARU, 
are just now coming on-line. Indeed, Yanagisawa & Muraki 
(1997), have done a preliminary study and conclude that the 
SUBARU telescope should be able to monitor about 200,000 
stars each in Pal 8 and NGC6453 with 10 second exposures. 
Longer exposures and stellar templates provided by HST 
may be able to push this number up considerably. Even if 
such large telescopes and high resolution are not available, 
application of the pixel lensing technique may allow monitor- 
ing of sufficient stars. Rhoads & Malhotra (1997) calculate 



that pixel lensing can allow 200,000 star equivalents to be 
monitored in M15 and an average over all globular clusters 
of about 60,000. We note that our clusters are considerably 
more luminous than the average cluster. Finally, it is impor- 
tant to realize that such experiments have already been per- 
formed! Mighell et al. (1992) reported two nights of obser- 
vations toward Omega Centauri with the Anglo- Australian 
telescope. Although no events were detected, 50,000 stars in 
a 12.5' field were monitored with 10 sec exposures. In an 
ongoing effort to search for variable stars, the OGLE col- 
laboration h as been examining po rtions of Omega Centauri 
and 47 Tuc ( Kaluzny et al. 1996 ). With 500 sec exposures 
in 1.6" seeing they follow about 30,000 stars in Omega Cen- 
tauri. In light of these numbers, the target figure of 100,000 
stars per cluster seems reasonable. 

The challenges involved in cluster microlensing searches 
should not be underestimated: following millions of very 
faint stars in extremely crowded fields scattered over a large 
part of the sky will require dedication, both of telescope 
time and its practitioners. The telescopes and techniques are 
both at the cutting edge of astronomy. However, optimism is 
warranted: 10 years ago microlensing was a pipe dream, now 
seven collaborations are actively pursuing it and millions of 
stars are monitored nightly. We must be bold in proposing 
and following up new ideas. 
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